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We show that integrable vertex and RSOS models with trigonometric Boltzmann 
Oh" 

weights and appropriate inhomogeneities provide a convenient lattice regularization for 



massive field theories and for the recently studied massless field theories that interpolate 
between two non trivial conformal field theories. Massive and massless S matrices are 
computed from the lattice Bethe ansatz. 

1. Introduction 

The perturbation of unitary minimal models Ait of conformal field theory with central 
charge 

6 , s 

c=l-- -, 1.1 

t(t-iy K ' 



by the (^13 operator of conformal weight 



fc = L r. (L2) 



Packard Fellow 



is known to give an integrable 1+1 quantum field theory ffl: A4A t . The infra red (IR) 



properties of this theory depend on the sign of the coupling. Writing the hamiltonian as 



the case g < gives the massive integrable field theory M.A t with a trivial IR fixed point 
while the case g > is "massless" and gives A4A^~ with an IR fixed point that coincides 



The conformal field theory Ai t is well known to describe the continuum limit of the 
(t — 2) th critical Ising model [[§, a critical point where in the Landau Ginzburg picture 
t — 2 minima of the potential coalesce. The massive flow describes the flow towards a 
phase where these minima separate and are degenerate and where spontaneous symmetry 
breaking occurs, while the massless one corresponds to one minima decoupling so in the 
IR the model is a multicritical Ising model with one order of criticality less. 

It is known that minimal models Ait of conformal field theory can be obtained by 
twisting of boundary conditions and then reduction of the Gaussian model with c = 1 
0. In a similar way, the massive regime of the perturbed theory with action (|1.3| ) can 
be described as a reduction of the sine-Gordon model || at coupling ^ = Due to 

twisting, the coupling constant in front of the cosine term scales as G oc (— g) 1 / 2 . The 
soliton mass m scales as m oc {—g)t. 

Formally the same correspondence between the two theories should also hold in the 
massless regime. In that case however the sine-Gordon coupling constant G turns out 
to be purely imaginary due to the above scaling. This is a regime that has not been 
much studied so far, and which presents subtleties 0. Notice of course that the rather 
violent non unitarity of this sine-Gordon model does not prevent the subset of observables 
corresponding to the minimal model to be unitary. 

Perturbations of coset models su(2) s ® sw(2)i_ s _ 2 /sw(2)i_ 2 (s = 1 above) by the 
operator of dimension h = ^p 2 - also give integrable models. We will denote them A4Af t 
depending on the sign of the coupling constant in front of the perturbation. These models 
can also be obtained by reduction of some other integrable models: the higher spin or 
fractional supersymmetric generalizations of the sine-Gordon model [[7] 

Further generalizations to coset models of type Gk X Gi/Gk+i are also possible. For 
each such model there is a natural integrable direction of deformation with the per- 
turbing field of dimension 1 — k J} +h v (where h v is the dual Coxeter number) leading to 




(1.3) 



withM t -i|l §• 



two different theories depending on the sign of the coupling constant. The role of the 
sine-Gordon model is then played by its Toda generalizations and further fractional super- 
symmetric generalizations [ fTU| . 

We show in this paper how inhomogeneous vertex or RSOS models with trigonomet- 
ric Boltzmann weights and appropriate inhomogeneities provide, in their scaling limit, a 
general lattice regularization of the above massive and massless field theories. This allows 
us in particular to put the massless scattering of |l l[] on firmer grounds, and to prove a 
conjecture in . In a way our approach is parallel to the one developed on [|12| for non 
linear sigma models. It is interesting that the lattice version of the sine-Gordon model 
which was studied in [O] turns out to be of the same nature. It was shown in |TJ[] that 
it can also be formulated as a model on the inhomogeneous lattice where at each vertex 
there is a cyclic representation of U q sl{2). 

Although the method is general, we discuss mainly the su(2) case. Appropriate ex- 
tensions to other groups are discussed in the last section. In section 2 we recall necessary 
results about the Bethe ansatz for inhomogeneous vertex and RSOS models. In section 3 
we discuss the thermodynamic limit. In section 4 we consider the scaling limit. In section 
5 we study the S matrix in the massless case and prove the conjecture of || . In section 6 
we discuss some examples of the method applied to other groups. 

Lattice regularizations of masssive field theories which are based on inhomogeneous 
integrable models with trigonometric Boltzmann weights have some advantages compared 
with the elliptic ones ||15|| . One is that in this case the lattice and continuum theory have 
the same symmetry. Another is that we know far more about the spectrum of trigonometric 
transfer matrices than about the elliptic ones. Finally for Lie algebras other than su(n) 
we do not expect elliptic vertex models JTB . 



2. Inhomogeneous lattice models of XX Z type 

2.1. Vertex models 

Consider a one dimensional lattice with N vertices formed by the intersection of N 
columns and one row. We assume the vertices are numbered from left to right. A U q sl{2) 
or generalized 6- vertex model on such lattice can be described as follows. Fix spins Sj, 
j = 1, . . . , N (in units where the fundamental has spin one) and an additional spin s. With 
each vertical edge we associate the space of states C Sj +1 and a complex number Uj. With 
each horizontal edge we associate C s+1 and a complex number u. 



Ui . . . Uj 



u 



Si ■ 



With each vertex we associate a system of Boltzmann weights which are given by R- 
matrices acting in the tensor product of corresponding spaces. In each space C s+1 we 
enumerate the basis by a = 1 , . . . , s + 1 . Fix an element of the basis in each space of 
states corresponding to the boundary edges, say a, a%, . . . , ajy, b, 61, . . . , 6jv- The partition 
function of the U q sl{2) model with the parameters and boundary conditions just described 
has the following form 



T s {u\{uMsj})lSi::t 



N 
N 



R ss i(u- Ul r c x---R ssN (u- UN )i* 



cjv-i«iv 



(2.1) 



ci,...,cjv-i 



Here R SSj is the U q sl{2) R-matrix described in JTTj and [fTS 

From this partition function with open boundary conditions on one row we can con- 
struct corresponding partition functions on rectangular lattices with various boundary 
conditions as traces of products of (|2.1| ) . 

Since the R-matrices satisfy the Yang Baxter equation we have commutativity of 
traces 



t 8 (u\{ Uj }\Uj}),t s (vliuMsA) 



0. 



(2.2) 



Here the transfer matrices t s (tt|{ , u :; }|{sj}) act in the tensor product 



vertex 



,SN 



c 



Sl+1 



(js N +i an j ^^gy are given by traces of ( |2.1|) with respect to "horizontal spaces" 

a=l 



(2.3) 



where d a = exp[/t(a — 1 — s/2)] for some complex number k. These transfer matrices 
commute with the total third component s z of the spin. The space H™ rte s JV we may call 
vertical space. We can thus regard ( [2.2|) as a generating function for a family of quan- 
tum integrable systems. The commuting family can be diagonalized simultaneously. The 
eigenvectors of t s (u\{uj}\{sj}) are parametrized by solutions {a/c} to the Bethe equations 



j=N l=M 

Y[ exp [ip Sj (a k - 2iuj)] ] [ exp [i<j>(a k ~ a t )} 

3=1 



■1, 



(2.4) 



1=1 



where M = • Sj — s z )/2 is an integer. The functions p s and <f> are given respectively by 



sinh i (a — si"f) 
exptp s {a) = — -f- — ■ — — , (2.5) 
sinh | (a + si^) 

and 

sinh^(a + 2z 7 ) 

expz0(a) = , f — - . (2.6) 

sinn ^{a — Ztj) 

With these conventions the dependence on heterogeneities in the vertical spaces appears 
only in the Bethe equations. Eigenvalues depend on the horizontal parameters u, s and 
the numbers a k as in the homogeneous case. Complete expressions, which involve a sum 
of s + 1 terms, can be found in [^9[. The dominant term in the thermodynamic limit reads 
(with proper normalization) 

Ajv'(«) « Y[exp{-i\pA<x k + 2iu) -tt]}, -> oo. (2.7) 



For more details and references see for instance []T9| . 
2.2. RSOS models 

If q is a root of unity, q = exp(i7r/t) where t is an integer, we can also construct the 
Restricted Solid On Solid (RSOS) models corresponding to the quantum group U q sl(2). 
States in such a model are associated with faces of the lattice, and they take values lj = 
l,...,t — 1. The adjacency conditions between neighbouring faces depend on the spin 
associated with the edge between them [^0[[^T1. The Boltzmann weights are associated 
with vertices of the lattice and they depend on the states of the neighbouring four faces. 



These weights are related to the vertex R- matrices via 6j calculus |22|||23|. 



Fix spins si, . . . , sn- The corresponding vertical space of states in the RSOS model, 

OS 
si,...,sjv 



Ti.ff os s with periodic boundary conditions, is best described by its basis 



{h, . . . , l N }, (k - k+i + Si)/2 E {0, . . . , sj, Si < k + <2t- Si. (2.1 



According to [21| , [24] the eigenvectors of the RSOS transfer matrix are obtained 



from the Bethe equations 



j=N l=M 

Y[ exp [ip Sj (a k - 2iuj)] exp [i(p(a k - a t )] = -w 2 , (2.9) 

3=1 1=1 



where we now have the constraint M = ^ ■ Sj/2 and u f = —(—) M y (Imw 7^ 0) and y is 
the eigenvalue of the operator given by 



7' 7' 

l\ ,...,ijv 



i=i 



that commutes with the RSOS transfer matrix. Eigenvalues of the RSOS transfer matrix 
are given by a sum of s + 1 terms as in the vertex case, each term being however weighed 
by different powers of u. For complete expressions see [E3| . These eigenvalues are a subset 
of the eigenvalues of the vertex model with u = exp/t (2TB) and in the spin s z = sector 
121 



2.3. Arrays 

We shall consider further only two choices of parameters s\, . . . , sjv and u\, . . . , u^: 
Case I: we choose si = . . . = sn = s and — u^j = «2,j-i = «A/2 for j = 1, . . . , N/2. 

Case II: we choose s^j-i = «i, S2j = S2 and — u^j = 1*2.7-1 = iA/2 for j = 1, . . . , iV/2. 

In both cases we assume N to be even and we consider transfer matrices as generating 
functions for hamiltonians of quantum spin models (of generalized XXZ or RSOS type). 
For hamiltonians we choose 



Hj = --4~ In \t s (iA/2 + u)(t s (-iA/2 - w))" 1 ! 
t au 



u=0 



and 



H n = -- 4- In \t s UiA/2 + u)(t s *(-iA/2 - u))- 1 ] 

-f- L ' J 



(2.10) 



(2.11) 



u=0 



t du 

so we have a total of four possible cases. 

In all cases the models consist of an even and odd sublattice and hamiltonians commute 
with the total translation operator T : j — > j + 2. Hamiltonians are local and can be written 



H = H+ + H 



(2.12) 



with 



H+(A) 



+2,j+i(*A)) 1 ^ +a ,j(0)i2 J - +2 j+i(iA) 

l_j odd i even 



(2.13) 



and 



#_(A) = -- Rj-iA-^)+ E ^i-2j-i(-iA)i?i- 2 j(0)(i? i _ 2J -i(-iA)) 1 

|_j even ^ Q dd 

(2.14) 

where dots stand for and R-matrices implicitely carry spin upper indices in agreement 
with the arrays I and II, and correspond to XXZ type or RSOS models. 

In terms of Bethe ansatz solutions eigenenergies of these hamiltonians read from ( |2.7| ) 

2 x - 

E I = -j2_^p s (a k + A)+p s (a k -A) (2.15) 

k 

and 

2 

E H = --J2 P S1 K + A) + p S2 (a k - A) (2.16) 

1 k 

(this corresponds to antiferromagnetic interactions. For the ferromagnetic case switch the 
sign of E). The momentum of the corresponding state reads 

Pi = J2p s (a k + A)+ Pa (a k -A) (2.17) 
k 

and 

Pi i = ^p Sl (a fc + A) +p S2 (a k - A) (2.18) 
k 

2.4- Physical comments 

Consider first array I. Call A the spatial period of the system (equal to two lattice 
spacings in the original square lattice). Call r the translation operator j ; — > j '< + 1 so r 2 = T. 
Then t s (iA/2)r and r(t s (— zA/2)) -1 describe vertical propagation as in the following figure 



while t s (iA/2) and (t s (-zA/2))" 1 describe propagation along the right and left diag- 
onal. The system can thus be intuitively considered as a "lattice light-cone" where 
the vertices represent interactions between s + 1 "bare" particles (see |26] for a re- 



lated approach). The natural hamiltonian associated with this description is given by 



Hi = i In [t s (iA/2)(t s (— zA/2)) _1 ] . It is however non local, and this causes some diffi- 
culties (see later). Array II would have a similar interpretation but with the left (right) 
diagonal carrying spin si (^2). Note that in this case the only interactions are between left 
and right moving bare particles. Let us emphasize that in such a picture we are dealing 
with quantum mechanics at real time. 

In the limit A — > 00 the first term in ( |2.13| ) and ( |2.14| ) decays exponentially while 
the second one remains of order one. In that limit, the second term in ( |2.13| ) and ( |2.14| ) 
describes a RSOS hamiltonian with some special boundary conditions, coupling respec- 
tively only even (odd) sites. In the limit A — > 00 the hamiltonian (|2.12| ) is thus the sum 
of two decoupled RSOS hamiltonians. When A is finite, they interact due to the first 
terms in ( |2.13| ) and ( |2.14| ) . It is easy to see on the other hand that the transfer matrix 
t s (iA/2)(t s (— iA/2)) -1 becomes trivial in the limit A — > 00. 



3. Thermodynamic limit 

We now consider the thermodynamic limit N — > 00. For the computation of thermo- 
dynamic quantitites, it is expected that the only relevant solutions of (|2.4p are given in 



terms of strings ||27|| . For technical simplicity we will restrict in the following to values of 



the anisotropy parameter q = exp(«7), where 7 = ir/t, t being an integer. 

3.1. XXZ type quantum chains 

For the vertex or generalized XXZ model case the solutions of ( |2.4| ) are made of n 
strings of length n = 1, . . . , t — 1 and of the (n = 1~) antistring of length 1 |27] (recall that 



an n string is a set of n complex numbers ojL = + ilj with / = n — 1, n — 3, — n + 1 
and an antistring of length one is a complex number with imaginary part equal to n). 
Regrouping the numbers a into their various strings leads to equations similar to ( |2.4|) but 
for the (real) centers of the strings 

j=N 

Y[ exp [ip n ^ s .{al - 2iuj)] JJexp [i(f) n , m (al - af 1 )} = -1, (3.1) 

j = l ml 

where the second product is taken over all possible strings and then over each center. 
In the thermodynamic limit equation ( |3.1| ) gives the following integral equations for the 
densities p n , p n of strings and holes 

5^Pn, Sj *V aj +^2<fin,m* Pm = 27Te(/9 n + p n ). (3.2) 



where V s encodes the distribution of the heterogeneous spectral parameters for the spins s 
(eg a sum of two delta functions or a single delta function in the array I or II) and e = — 1 
for the antistring. 

In the following we set a n>s = p n ^ s and A n m = 2ir8 n>m 5 — </> n>m . The various above 
functions are more conveniently expressed in terms of the Fourier transforms of their 
derivatives. For n,s= 1, . . . , £ — 1 one has 

sinh(nx) sinh[(£ — s)x] 
ttn,s — . r 77 \ ■ i / \ j n < s. (3-3) 
smh(£xj sinh(x) 

and 

? cosh(x) sinh(nx) sinh[(£ — m)x] ,„ 

A n , m = 2 . r 77 r ■ w \ — , " < ™- (3.4) 

smh(£xj sinh(x) 

when one or two of the labels involve the antiparticle one has 

sinh(sx) ,„ r\ 

a l~,n — • 7" 77 \) 

smh(£xj 

and 

i , , . sinh(nx) . „ . .„ „. 

smh(£x) 

~ sinh[(£ — 2)xl - * . „. 

Ai-.t-i = . L Y ; J , A!- fl - = A ±1 (3.7) 

smh(£x) 

and the kernels are symmmetric otherwise. 

Once the thermodynamic limit has been taken we compute the density of free energy of 
the infinite chain (with still finite lattice spacing) T = £ — TS, where T is the temperature 
and S, S are the densities of energy and entropy respectively. The formula for the entropy 
is the iV — > oo limit of the combinatorial entropy of distribution of holes and particles |[28|| 



/oo 
[(Pn + Pn) log(Pn + Pn) ~ Pn log(p„) - pn log(/5 n )] • (3.8) 
-°° 

The expression for the energy depends on the particular array. For array I one has 

2 f°° 

Si = -- / [o„, s (a + A) + a n , s (a - A)] p n {a)da, (3.9) 

and for array II 

2^ r 

£11 = -- y I [a n , Sl (a + A) + a UjS2 (a - A)] p n (a)da, (3.10) 

" „ J — oo 



3.2. RSOS quantum chains 

We can follow the same steps for the RSOS chain. The additional phases in 
seem to disappear in the equations for densities in the thermodynamic limit, but in fact as 
discussed in , they suppress the t — 1 string and the antistring. This implies moreover 
that for the t — 2 string there are no holes. One can then eliminate pt-% in the various 
equations to get a system exactly identical to ( |3.2| ) but where n = 1, . . . , t — 3 and where 
t is replaced by t — 2 in the various kernels 



3. 3. Thermodynamic equilibrium for XXZ type quantum chains 

The minimization of the free energy leads to coupled equations for the densities |27| 



|2l| . For array I we have the following. Introducing e n = Tlog(p n / p n ) we get, after 
inversion of the matrix in ( |3.2| ) , 

-e(a)* nia = e n -ra*^JV^log[l + exp(e m )/r]. (3.11) 

m 

where we have defined 

e(a) = s(a + A) + s(a - A). (3.12) 

and 

s(a) = — — —-. (3.13) 

v ; cosh(ta/2) v ; 

In the last expression N** z is the adjacency matrix of the following diagram 

O t-1 
12 s t-3 I 

0--O -O- -0--0 t-i 

For the free energy per unit length we find 

/oo 
{Tlog[l + exp(e s /T)] + a s , s (a)}e(a)da. (3.14) 
-oo 

Similarly for the array II, inversion of the matrix in ( |3.2j ) now leads to 

-e Sl (a)5 nySl - e S2 (a)5 njS2 = e n - Ts * ^ N n,* Z ^gl 1 + exp(e m /T)], (3.15) 

m 

where 

e Sl (a) = s(a + A), e S2 (a) = s(a — A), (3.16) 

and for the free energy 

F n = - J {Tlog[l+exp(e Sl /T)]+a SlS Je Sl (a)da+{Tlog[l+exp(e s ^T)]+a S2S2 }e S2 (a)da. 

(3.17) 



3.4- Thermodynamic equilibrium for RS OS type quantum chains 



Things are very similar in that case. The only difference is that sums over solutions 
run only over n = 1, . . . , t — 3 strings in ( |3.11| ) , ( |3.15| ) , (|3.14| ) and ( |3.17| ) while in ( |3.11| ) 
and ( |3.15| ) N^^ s is the incidence matrix of the diagram 



o- 



2 

o 



-o 



t-3 

-o 



3. 5 T — ► limit for array I 

Consider first the ground state of array I at T = (in this paragraph the distinction 
between generalized XXZ and RSOS chains is simply encoded in the two possible incidence 
diagrams). System ( |3.11| ) reads then 

-e(a)5 ntS = e n - s*2J N n, rn et l , (3.18) 

m 

where e + (x) = e(x)H(e(x)) is the positive part of the e function (H being the Heaviside 
step function). The solution of this system is e+ = and e~ = — ed ns . From (|3.18|) we see 



that the ground state is filled with s strings. Excitations above this ground state are made 
of holes in the sea of s strings, with excitation energy precisely equal to e (other strings 
are important in determining for instance the spin of the excitations, but they do not 
contribute to the energy and momentum). The momentum of the excitations is similarly 
found to be 

p = 2 tan" 1 ( Sm ^ ta / 2 ) (3.19) 



and as in [[29| they always occur in pairs (with our convention p is defined modulo 4ir). 
There are two Brillouin zones. In the first one the model has massive excitations with a gap 
at p — 77 that vanishes exponentially with A. In the second zone the model has massless 
excitations around p = 0, 2n. This is true both for XXZ and RSOS chains provided 
s < t — 3. In our subsequent discussion we discard the finite parts of the momenta. The 
role of these contributions in the scaling limit is similar to the one in spin s = 1 case 
discussed eg in [[30] 



3.5. T — > limit for array II 

In that case the ground state is filled with si and S2 strings and the excitations above 
it are made of holes in the sea, with excitation energies equal to e Sl , e S2 . The momentum 
of excitations is similarly found to be 

p Sl = 2tan _1 {exp[t(a + A)/2]}, p S2 = 2 tan _1 {exp[t(a - A)/2]}. (3.20) 

There are massless excitations around p = 0, n in both cases. Again this is true for both 
XXZ and RSOS chains provided si, s% < t — 3. 

4. The scaling limit 

4-1. XXZ type chain and array I 

In the first Brillouin zone the dispersion relation exhibits a gap that vanishes in the 
limit A — > oo as 

mA = 5E w 4exp(-tA/2), (4.1) 

Let us consider the scaling limit. Let >> 1 and at the same time A << 1 in such a 
way that the dimensional length of the chain L = AN is much bigger than the correlation 
length ^ = jg, the latter itself remaining finite. Hence we need 

A,iV >> 1, A « 1, NSE » 1, A oc 5E (4.2) 

We then concentrate on the contributions of order A to the hamiltonian. At finite rapidity, 
that is very close to the origin of the first Brillouin zone, we then find for the excitation 
energies 

e(a) w 4e~ tA/2 cosh(to/2) = Am cosh 9, (4.3) 
where the mass m = l/£, and the physical rapidity 6 = ta/2. Similarly for the momentum 

p{ot) w 4e" tA/2 sinh(ta /2) = Am sinh 6, (4.4) 

Relations ( |4.3| ) and ( f4.4j ) correspond to relativistic particles of mass m. Excitations at finite 
distance from the origin in the first Brillouin zone have finite energy and will disappear in 
the scaling limit. 



In the second Brillouin zone we also have excitations with energy of order A. As 
a » A (A being itself very large) we have 



e L la) 



2e tA/2 e -ta/2 



(4.5) 



with 



Similarly as a « —A we have 



Pl{oc) « -e L . 



e R (a) « 2e tK ' 2 e ta '\ 



(4.6) 



(4.7) 



with 



p R {a) « ei?. 



(4.f 



As in |29| there are only even number of the different types of excitations and the finite 
contributions to the momenta can be discarded. These relation correspond to massless left 
(respectively right) moving particles. 
As A — ► we expect the following 



H, 



latt 



N£ + H scal A, 



(4.9) 



where H sca i is the hamiltonian of the continuum field theory obtained in the scaling limit. 
On the other hand from |JT[], we observe that the thermodynamic equations ( ^.11| ) of 
our system, for energies of order A, coincide in the first Brillouin zone with the so called 
Thermodynamic Bethe Ansatz equation J for the spin s or fractional supersymmetric sine- 
Gordon model with coupling constant 



2 



SG 



8n s(s + t) 



(4.10) 



where t = t — s and cjjv = ^pf? once the temperature in our problem has been identified 



as usual with the inverse radius R in , [[32] . Similarly in the second Brillouin zone the 
thermodynamic equations of our system for energies of order A coincide with the " massless 
TBA" equations for the generalized Gaussian model (a system made of a free boson and 
Z s parafermions) at radius of compact ificat ion appropriately related to (|4.10|) |33|]. The 



2 Here by TBA we refer to thermodynamic analysis of a gas of relativistic particles in one 
dimension, with factorized scattering. 



excitations are then massless since one is describing a conformal field theory. Moreover, as 
A >> 1, the regions of interest in the first and second Brillouin zone completely decouple 
since the kernel in the lattice equation (|3.11| ) decays exponentially outside a region of finite 



size in a variable. We therefore conclude that 

H scal = H SG -s + H G _ s , (4.11) 

where the first term is the hamiltonian of the spin-s sine-Gordon model and the second 
term the hamiltonian of the generalized Gaussian model (a conformal invariant theory) on a 
cylinder of circumference L with appropriate boundary conditions: Hq = ^-(Lq + Lq — jo)- 
The identification of the scaling limit could be made completely by computing, from 
the lattice model [19| [ 34 1 fl35 1 , the S matrix for the scattering of excitations. For doing 



so one interprets the holes in the sea as the physical particles. The other string solutions 
correspond then to pseudoparticles. In the present case three S matrices would have to 
be computed: the one describing the massive SG model and two others, identical ones, 
describing scattering between left (respectively right) "massless particles" in the conformal 
field theory. One would thus get three TBA diagrams like 



12 s 

o- -o- -®- 



where the marked node corresponds to a massive particle of mass m , repsectively to a left 
or right moving particle. We comment more on massless scattering in section 5 where we 
also give a detailed example of an S matrix computation. Let us just mention that the S 
matrices extracted in this fashion coincide with the known one for spin s sine-gordon 
and the ones for generalized Gaussian models [[33] |29| [j36 |. 




To understand more directly what perturbed conformal field theory we are dealing 
with we finally consider [|37|] a limiting process similar to (fO|) but with the correlation 



length much bigger than L ie 

A, N -> oo, A -> 0, N5E « 1, A oc 5E (4.12) 



In that limit we have 

Hscal = Hconf + ©Xp(— A)H pert , 



(4.13) 



where H con f = Hq_ s + Hq- s (this corresponds to ( |2.12| ) being a sum of two decoupled 
XXZ type models in the limit A — > oo) and H pert is determined by the first term in ( 2.13Q 



and ( 2.14Q . Comparing this with the formula giving the physical mass leads to the value 



of the conformal weight of the perturbation h$G-s — This conformal weight could 

also be found for s = 1 by noticing that taking the limit A — > oo the equations ( |2.4|) with 
s = 1 reduce to the ones of the Thirring model with bare mass uiq = 4sin7exp(— A) and 
bare rapidity a |34[] [[381 . Notice that the conformal weight does not depend on the spin 



s, ie it depends only on the skeleton of the TBA diagram and not on which node(s) are 
massive. This result can also be directly established at the level of the TBA equations 



11 . 



4-2. RSOS chain and array I 



Things are quite similar in that case, with the additional truncation of (|3.11| ) . The 
thermodynamic equations (|3.11|) for energies of order A now coincide in the first Brillouin 
zone with the TBA equations for the theory A4A~ t written in || . Similarly in the 
second Brillouin zone the equations coincide with the massless TBA that describes minimal 
conformal field theories |3^| . These three TBA are associated with the same diagram 

12 s t-3 

o — o <g) -o 

where the marked node corresponds to a massive particle of mass m, respectively a left or 
right moving particle. Equation ( 4.11| ) now reads 



Hscai = H MA - + H Mst , (4-14) 



where the notations are explained in the introduction. In the limit ( |4.12| ) one finds 

H sca i = H conf + exp(-2A)H pert , (4-15) 

where H con f = Hj^ st + Hj^ st , H pert is determined by the first term in ( |2.13|) ( |2.14| ) . 
The main difference with the unrestricted case discussed above is that the perturbation 
starts with order exp(— 2A) instead of order exp(— A): this is due to the different spectral 
parameter dependence of the R-matrices (the different gradations from U q sl(2) point of 
view 1 13 911)1. As a consequence one finds h = as announced in the introduction. 



t 

To conclude this subsection we have established that array I provides a regularization 
of the interpolating theory A4A~ t (plus a "background" M. s t theory that decouples). 



The difference of exponent in ( |4.13| ) and ( 4.15 ) can also be observed in the ellpitic case: see 



the Boltzmann weights in [21| . 



4-3. XXZ type chain and array II 

We consider again the scaling limit defined in ( |4.2| ) . In the region \a\ >> A (the 
second Brillouin zone) things are the same as for array I. They are quite different in the 
first Brillouin zone. In the present case we have, in the limit A — > oo, and for finite 
rapidities, 

Am a Am a 

where m, 9 are the same as before. It is straightforward to study as well the momentum 
of the excitations and one finds 

p Sl = -e Sl , p S2 = e S2 , (4.17) 

ie the dispersion relations for left and right moving massless particles. We thus have left 
and right moving massless particles in both Brillouin zones. However, while in the second 
zone there is no left-right scattering since the two types of particles are widely separated 
in the rapidity variable (a » A and a « —A), in the first zone we expect to have non 
trivial left right scattering. This leads to a theory that is non scale invariant. In general 
however, one does not get a meaningful scaling limit due to the lack of symmetry between 
left and right sectors. The first step to recover this symmetry is to consider the RSOS 
version. 

4-4- RSOS chain and array II 

In the scaling limit, thermodynamic equations in the first Brillouin zone coincide 
with then with the TBA equations associated with a diagram of type A n _3 with physical 
particles at nodes s± and S2- For LR symmetry one then has to require s± = s = t — 2 — 82, 
leading to the TBA diagram 

1... s... t-2-s... t-3 

O- - -<g) <8>- - -O 

m — 9 m 9 

2 2 

This TBA coincides then with the one conjectured in || , [f|D[ for the flow between 
successive minimal models. Hence we expect 



H SC al - H MA+ + H M st • 



(4.18) 



Observe on the other hand the following important symmetry for a RSOS model with 
t — 1 heights (with appropriate normalization) 



7T 



(this can be derived from identities in JT^] together with the observation that in the latter 
reference the spectral parameter u in the weights W pq is shifted, when compared to ours, 
by the amount ^t^t). 

Hence array II in the RSOS case is equivalent to array I, but with the cutoff A acquiring 
an imaginary part i^. We can now make the identification with perturbed minimal models 



more complete by considering the limit ( 4.1 2| ) . Due to the imaginary part of A we observe 



indeed the change of sign of the coupling constant in (^4.15| ) for the RSOS model: 

H sca i = H Mst + H Mst - exp(-2A)H pert , (4.20) 

with Hp er t the same as in ( f4.13| ) . These results indicate that we have a regularization of 
the interpolating theory A4A^ t . To complete the identification we shall compute the S 
matrix in the following section. 

Notice that, based on Q4.19f) , it is tempting to consider also a type I array for the 
vertex model with a cutoff A having imaginary part ^ . From ( f4.13| ) it corresponds indeed 
to the spin s sine-Gordon model with imaginary coupling constant. More detailed study 



of this model should be of interest [41 



4-5. Computation of central charges 

In the framework of the TBA it is standard to compute the central charge in the UV 
and IR. Let us notice that it is also possible to extract these two numbers directly from 
the lattice model by considering appropriate limits. In the scaling limit we expect the 
entropy <S(T, A) to become a function of the scaling variable — where r = ^. The UV 
limit corresponds to r >> m and the IR limit to r « m. To extract the central charge 
from the lattice model one has to proceed as follows 

rT 



and 



5(T, A) « ^-c IR if T « e" tA / 2 , T - 0, A - oo, (4.21) 
ttT 

S(T, A) w — cuv if T » e"* A / 2 , T -> 0, A -> oo, (4.22) 
3 

For the RSOS model and array I the results of this computation are cjjv = 2 ^1 — ^ryj j 
and cm = 1 — ^|z?n wnue f° r array II cuv is the same but cj^ = 1 — ^z?n + 1 — (t-i)(t-2) • 



5. The S matrix for massless flows between minimal models 

5.1. Massless S matrices 

To start we briefly discuss the definition of massles S matrices. For integrable sys- 
tems one can define scattering amplitudes via symmetry properties of eigenstates of the 
hamiltonian of the system under permutation of particles. We parametrize energy and 
momentum of left and right moving massless particles as in Q4.16|) and fl4.17|) so if ^qr ql 
is an eigenstate of H and P one has 

and 



* J 
J-R t)R qL nL 



Consider a given set of rapidities 9 H , . . . , 9 N , 6>f , . . . , 9 M and order them O^m < . . . < 
9 a (N+M) where a E Sn+m- To a given order characterized by some element r of Sn+m 
one can associate a linear space H. T (9 R , 9 L ) of dimension C d k n + m ) where D is the number 
of isotopic degrees of freedom. All these spaces are isomorphic and the isomorphism 

H T (9 R ,9 L )*H T ,(9 R ,9 L ), r,r' e S n +m (5.3) 

defines some S matrix S TjT >(9 R ,9 L ). The two particle S matrix is obtained in the case 
t = a and r' = ooi where Oi transposes i and i + 1 as 

Server = Sr R (9 r - 9 R ) if CT (j) = 9 R , CT (i + i) = 9 R , (5.4) 

S a ,ai(j = Sll(9j — 9%) if 0cr(i) = 9 j , = 0jt , 

and 

Serbia = Slr(9^ — 6*^) if = 0^, 0(r(i+l) = 0jf , 

with a similar equation for S^z,. 

When we put the theory on a line of length L this symmetry implies the form of the 
eigenfunctions 

^0r >9 l(x 1 ,...,x n+ m) = ^rexp[z(p <T ( 1) xi + . . . +P<t(n+m)Xn+m)], (5-5) 



where the amplitudes A T behave as ( |5.3| ) . The standard form of periodic boundary 
conditions |33| then follows. 



5.2. S matrix for M.A^ t 



To get the S matrix of the excitations we proceed as in |35[] . For simplicity we 
first consider the case s = 1. Start from 

^ t-3 
[ a^ 2 \a + A) + a%-%(a - A)l = 2np n + ^ * p m , (5.6) 



2 



n 



m=l 



where the subscript mean £ is replaced by £ — 2 in equations ( |3.3| ) up to ( |3.7| ) . Rename 
now <jl = pi, <jr = pt-3. The densities o~l,r now are considered as describing physical 
particles (they are densities of holes in the original Bethe ansatz due to the structure of 
the ground state and its excitations). Similarly rename a n = p n -x, = ex, en = et_3 
and isolate the physical densities. After lengthy calculation one finds 

m=t— 5 

27ra n + A^; n 4) * cr m = c4 t_4) * a L + a^5 5 * <r R (5.7) 

as well as 

n=t— 5 

2yr(a L + o- L ) = e L + af~ 4) * a L + af_" 5 4) * a R - ^2 a n~ 4) * a n, (5.8) 



m=l 



n=l 



and 

n=t — 5 

27r(cr fl + Oi?) = e fl + af_~ 5 4) ★ a L + aj i_4) ★ a fl - ^ ojt^i * cr n , (5.9) 

n=l 

here a, A refer to the same kernels as in (|3.4|) ( |3.7|) but with £ replaced by £ — 4. As a 
consequence there is L : R symmetry ie An,m^ = 4-2-n t-2-m 1 Taking now the limit 
A — ► oo, the equations (|5.7| ) ( |5.8| ) ( |5.9| ) can be identified with the ones one would write for 
a theory with L particles of spin 1, R particles of spin £ — 3 with factorized scattering given 
by Sll = <S^3) Srr = •S'tZa'* -5 = Sll and = S]^ 5 @, where in each case S is made 
of Boltzmann weights of the appropriate RSOS model with £ — 3 heights! and adjacency 
conditions depending on the spin, and the usual minimal prefactors ensuring unitarity and 
sing symmetry i. Indeed (|5.7|) just expresses the densities of pseudoparticles in the 



crossing 



problem of computing the eigenvalues for passing a L or R particle through a set of L and 



4 These identifications hold up to a phase, which cannot be determined from the TBA equations 

5 To be put in correspondence with the minima in the Landau Ginzburg picture 

6 These identifications hold up to a phase, which cannot be determined from the TBA equations 



R particles with the densities cr^, cxr (that is the eigenvalues of the "monodromy matrix" 
T s in section 2). Similarly ( |5.8|) and (|5.9|) express the periodicity condition for the wave 
function of L and R particles. 

Using the symmetry between spin 1 and spin t — 5 in the RSOS version, one can 



rewrite Slr(9) as S 



ii 



as well, in agreement with the result conjectured in 
The computation easily extends to the higher s case. One finds then 

S LL = Sl 1 ^ ® s t-i-2s ® 1, Srr = 1 ® S t -i-2, ® (5-10) 

and 



S M = 1 ® ^i_ 2s 



<g>l, (5.11) 



2 

where S n is the S matrix based on the RSOS model with n heights. In the above notation, 
left and right components of the tensor product correspond respectively to L and R RSOS 
models (with s + 1 heights) , while the middle component is as in the case s = 1 above, ie 
it contains both L and R RSOS models (with t — 1 — 2s heights), and the S matrix acts 
on LL, RR, LR respectively. 

5.3. Physical comments 

Notice that in the interacting theory we have non trivial LL,RR,LR scattering, while 
in the lattice light-cone picture only LR scattering occurs between the bare particles. 

It may be useful to summarize the various correspondence between lattice RSOS 
models and thermodynamic equations, RSOS models in S matrices, and TBA. In all the 
previous sections the lattice model is based on quantum parameter q = exp(z7r/t). This 
corresponds to a lattice RSOS model with t — 1 heights. In the massive case, the continum 
limit is characterized by a quantum parameter q m assive = exp(z7r/(t — 1)). The S matrix 
is based on a RSOS model with t — 2 heights. The TBA diagram has t — 3 heights ( t — 4 
come from pseudoparticles, one is massive). In the massless case the continuum limit has 
parameter q maS sless = exp(in/(t — 2)). The S matrix is based on a RSOS model with t — 3 
heights. The TBA diagram has still t — 3 heights (t — 5 come from pseudoparticles, two 
massless). For instance it we take t = 5, s = 1 that is the lattice A4 model, in the case of 
array I we get the massive flow from the tricritical Ising theory (c = 7/10) to the trivial 
fixed point, and in the case of array II the massless flow from tricritical Ising to Ising. 



6. Generalizations 



The method can be extended easily to more complicated arrays, or to other algebras. 
As an example let us focus on so(2r). As is generic in the higher rank case, the analysis of 
the Bethe ansatz equations, together with a string hypothesis for every color a = 1, . . . , r 
leads to a TBA system that can be associated with a diagram in the plane order of roots 
versus color [|37fl . Restricting to the case q = exp(z7r/2r) and turning moreover to the RSOS 



case there are only real solutions (1 strings) to take into account here, and the diagram 
is again a D diagram, where this time the nodes are labelled by colors. We can then 
consider several kinds of arrays as before. We restrict now to ferromagnetic interactions 
@. A simple case is the type I array built with the fundamental representation. Then one 
recovers the same results as in above (see discussion before eq (|4.10| ) for s = 1), due to 



the well known fact that the rational points of the sine Gordon model /3 2 = 8n^^ can 
also be described as cosets so(2r)i <g> so{2r)i/ so{2r) 2 (a similar correspondence with the 
sine-Gordon model at (3 2 = would be observed in the antiferromagnetic case, where all 
nodes would be massive, corresponding to the bound states). This correspondence extends 
to higher spin. The most interesting case occurs with the spinor representations. Consider 
first the array of type I based say on positive chirality. One gets then in the limit A — ► oo, 
thermodynamic equations that coincide with the TBA associated with 

m cosh 9 

1 r-3 j 

O- -O- -O- -O r - 2 



This corresponds to the massive perturbation of Z r parafermions by the operator of 
dimension h = ]43|]. For the array of type II, using the two chiralities for s± and s 2 one 



finds in the limit A — > oo thermodynamic equations that coincide with the TBA associated 
with 

1 r-3 

0--0- -o- O ' 




7 Recall that in level rank duality, ferro and antiferro interactions are exchanged [|l9| [42] 



which corresponds to the massless perturbation of parafermions with cuv = ^7x1 5 c iR = 
1~ ( r+ i) 6 ( r -|-2) • ^ we now consider the r — > 00 limit, these two possibilites provide integrable 



discretizations of the 0(3) non linear sigma model at G = (resp. G = it) following [fJ3 



Finally, it is easy to read the S matrix for these two problems. For doing that we 
interpret the TBA diagram again as an sit (2) diagram. In the first case the A r -i part 
then arises from the pseudoparticles and the additional leg indicates that the physical 
particles have spin 2. Hence S = S^L X . In the second case we have to interpret the A r -2 
part as related to pseudoparticles, and the two additional legs correspond to the physical 
L and R particles. Then Sll = Srr = Srl = Sy-2- 



7. Conclusion 



In conclusion we would like to notice that the approach where one chooses for hamil- 
tonian (eg in the type I array) 

Hi = \ In \t s (iA/2)(t s (-iA/2))- 1 ] 
1 L J 

(as is done in |2^] ) seems correct for the study of the scaling limit but breaks down in 
the perturbative analysis ( |4 . 1 2|) . In that case instead of equation ( |4.13j ) or ( |4.14| ) one gets 
a trivial unperturbed term, and the identification with perturbed conformal field theories 
becomes questionable. This problem is presumably due to the hamiltonian being non local. 
We also would like to comment about the study of higher conserved quantities from 



the lattice. For simplicity we restrict to the massive case. Expand the energy (|3.12| ) in the 
region |a| < A, after renaming e — > e^ 1 -*, 



n=0 



(2n + 1 



tk 



cosh 



(2n + 1 



ta 



(7.1) 



Instead of ( |2.1U| ) we can consider other hamiltonians obtained by taking higher derivatives. 
Define therefore 



H 



(n) 



du 7 



In [t s (iA/2 + u)(t s (-iA/2 - u))' 1 ] 



u=0 



Its eignvalues are given by 



(7.2) 



E i U) = ("^E^W A) - (-l) n #(«fc-A). (7.3) 



The analysis of thermodynamics and excitations is very similar to the case n 
finds the energies of excitations to be given by 



e (n) («) 



n-l 



d 
dX 



n-l 



,(1), 



a). 



1. One 



(7.4) 



Expand now these quantities as (|7.1[) . One finds for instance 



e (2) = 4^(2n + l)exp 



n=0 



J A' 
-(2n+l)y 



cosh 



(2n + l) y 



(7.5) 



All these quantities have dominant behaviour e™(a) ~ 4e tA / 2 cosh(ta/2). However by 
forming linear combinations of them we can extract higher order terms. For instance 



e (2)_ e (l)^ 8e -3tA/2 cosh(3to/2) 



(7.6) 



Hence 



(7.7) 



with obvious generalizations to higher orders. In the same way as H^} al is the quantum 
sine-Gordon hamiltonian, the H^J^ 1 ^ are the higher hamiltonians - non trivial for odd 
grade only - of the integrable hierarchy [|4 . They commute since the various latttice 
hamiltonians ( |7.2| ) used in their construction do. 
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